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Abstract 

We investigate numerically the Self Organized Criticality (SOC) properties of the 
dissipative Olami-Feder-Christensen model on small-world and scale-free networks. 
We find that the small-world OFC model exhibits self-organized criticality. Indeed, 
in this case we observe power law behavior of earthquakes size distribution with 
finite size scaling for the cut-off region. In the scale-free OFC model, instead, the 
strength of disorder hinders synchronization and does not allow to reach a critical 
state. 
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1 Introduction 

The idea of the seismogenic crust as a self-organized complex system was 
introduced over the years as a possible explanation for the widespread occur- 
rence of space-time long-range correlations in earthquakes dynamics, similar 
to those observed in critical phase transitions [1]. In general, the term self- 
organized criticality (SOC) [2] refers to the intrinsic tendency of a large class 
of spatially extended dynamical systems to spontaneously organize into a dy- 
namical critical state. One signature of SOC is the presence of both a power 
law behavior in earthquakes size distributions and a finite size scaling for their 
cutoffs. Among the great number of different SOC models [3,4] developed in 
the last years, the OFC model [5] , introduced by Olami, Feder and Christensen 
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in 1992, has played a key role in modelling earthquakes phenomenology. How- 
ever the presence of criticality in the non-conservative version of this model 
has been controversial since its introduction [6] and it is still debated [7,8], 
also in relation with the influence of topology. In literature, OFC models on 
different topologies has been investigated, in particular the 2D nearest neigh- 
bor lattice (NNL) model [9], annealed random neighbor (ARN) graph model 
and the OFC model on a quenched random (QR) graph [10]. The purpose 
of our work is to study the effects of small-world (SW) and scale-free (SF) 
topologies on the criticality of the non-conservative OFC model. 
The paper is organized in the following way. In section 2 we review the OFC 
model and we point out the main reasons that have induced us to study the 
non- conservative OFC model on SW and SF topologies. In section 3 we in- 
vestigate the SW OFC model: in subsection 3.1 we show the earthquakes size 
distributions for the non- conservative SW OFC model and in subsection 3.2 
we characterize the critical behavior of the model through the finite size scal- 
ing ansatz. Finally, in section 4 we investigate the OFC model on a scale-free 
network, obtained by preferential attachment procedure [11]. Conclusions are 
drawn in section 5. 



2 The Olami-Feder-Christensen model 

The Olami-Feder-Christensen (OFC) model [5] is defined on a discrete system 
of N sites (blocks or of fault elements) on a square lattice, each carrying a 
seismogenic force (see Fig. 1). Such a force is simulated by associating to each 
site a real variable Fj, which initially takes a random value in the interval 
(0, F t h). All the forces are increased simultaneously and uniformly (mimicking 
a uniform tectonic loading), until one of them reaches the threshold value F t h 
and becomes unstable (Fj > F t h)- The uniform driving is then stopped and 
an "earthquake" (or avalanche) starts: 



where "nn" denotes the set of nearest-neighbor sites of i. The parameter a 
controls the level of conservation of the dynamics and, in the case of a graph 
with fixed connectivity q, it takes values between and 1/q (a = 1/q corre- 
sponding to the conservative case). The toppling rule (1) can possibly create 
new unstable sites, producing a chain reaction. All sites that are above thresh- 
old at a given time step in the avalanche relax simultaneously according to 
(1) and the earthquake is over when there are no more unstable sites in the 
system (Fj < F th , Vi). The uniform growth then starts again. The number of 
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topplings during an earthquake defines its size, s, and we will be interested in 
the probability distribution P^(s). In the following the boundary conditions 
of the model will be "open" , implying that F = on the boundary sites. 
At this point it is important to emphasize that the OFC model behavior 
strongly depends on the chosen topology. For instance, in the dissipative NNL 
OFC model with open boundary conditions the earthquakes size distribu- 
tion is described by a power law [9,12], characterized by a universal exponent 
t ~ 1.8 independent of the dissipation parameter. However, at variance with 
the conservative case where a full SOC behavior is observed, finite size scaling 
appears to be violated in the pdf cutoff-region (see Table 1). 




Fig. 1. Critical stress field of a 64x64 lattice (NNL OFC model) in the critical state. 

In ARN OFC models [13,14,15,16], where each site interacts with randomly 
chosen sites instead of its nearest neighbors on the lattice, there is criticality 
only in the conservative case, where it becomes equivalent to a critical branch- 
ing process. As soon as some dissipation is introduced, the earthquakes become 
localized although the mean earthquakes size diverges exponentially as dissipa- 
tion tends to zero and there is no power law distribution (see Table 1). Actually 
it is interesting to point out that criticality in the OFC model on a lattice has 
been ascribed to a mechanism of partial synchronization [17]. In general the 
system shows a tendency to self-organize into a periodic state [17,18,19] which 
is frustrated by the presence of inhomogeneities such as the boundaries. In 
addition, inhomogeneities induce partial synchronization of the elements of 
the system building up long range spatial correlations and a critical state is 
obtained. The mechanism of synchronization requires an underlying spatial 
structure and therefore cannot operate in an ARN model, where each site is 
assigned new random neighbors at each update. 

In the OFC model on a QR graph, where the choice of neighbors is not an- 
nealed but quenched and all the sites have exactly the same number of nearest 
neighbors q (both for q = 4 and q = 6), the dynamics organizes into a subcriti- 
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cal state. This is analogous to what happens in the OFC model on a NN lattice 
with periodic boundary conditions, where no critical behavior is observed at 
all [17,18,19]. In the QR case, in order to observe scaling in the earthquakes 
distribution, one has to introduce some inhomogeneities [10]. In particular, it 
has been found that it is enough to consider just two sites in the system with 
coordination q — 1 [10]. When either of these sites topple according to rule 
(1), an extra amount aFi is simply lost by the system. In such a way spatial 
correlations can develop, because the topology is quenched, there is power law 
in earthquakes size distribution and also finite size scaling is observed (see 
Table 1). 

In this work we study the non-conservative OFC model on both a small-world 
and a scale-free topology. 

First of all we expect that the inclusion of some inhomogeneities in the sites 
degree is not the unique way to obtain SOC. Indeed, as we are going to show, 
an alternative way is to keep fixed the sites degree and to change the topology 
of the underlying network, for instance by considering a small-world graph ob- 
tained by randomizing a fraction p of the links of the regular NN lattice. Here 
we will use the term "small-world" to refer to a rewired lattice (with fixed 
connectivity) with the minimum number of rewired links such that the char- 
acteristic path length L is almost as small as that one for the corresponding 
random graph [20,21]. As shown in the right panel of Fig. 2, this is obtained 
already for very small values of p (p ~ 0.01), much before the random graph 
limit (p — 1). 

A small-world topology is expected to be a more accurate description of a real 
system according to the most recent geophysical observations that indicate 
that earthquakes correlation might extend to the long range in both time and 
space [22,23]. In fact, if a main fracture episode occurs, it may induce slow 
strain redistribution through the earth crust, thus triggering long-range as 
well as short-range seismic effects [24,25,26,27,28]. The presence of a certain 
percentage of long-range connections in the network takes into account the 
possibility that an earthquake can trigger other earthquakes not only locally 
but also at long distances. 

On the other hand, one can consider a different kind of networks with a small 
L, the so called "scale-free" networks, which differ from the small-world graphs 
for having a power law distribution of the site degree. Scale-free networks are 
very common in nature and have also been used for SOC models (see ref. [29] 
for sandpile dynamics on SF network) but they have not been investigated, 
as far as we know, in the context of OFC models. It is known that, when the 
connectivity is not fixed at all but only in average (as for a random graph 
in Ref. [10]), the strength of disorder is enough to destroy critical behavior. 
Thus we expect that for SF networks, where the connectivity has a power law 
distribution, synchronization will not take place and it will be not possible for 
the OFC model to reach a critical state. In the last part of the paper we will 
show that this is exactly what happens. 
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3 The OFC model on small-world network 

To investigate the effects of the small-world topology on the criticality of 
the non-conservative OFC model, we follow the method proposed by Watts 
and Strogatz to construct a network which interpolates between a square NN 
lattice and a quenched random graph [20,30]. We start with a two-dimensional 
NN square lattice in which each site is connected to its 4 nearest neighbors. 
The links of the lattice are rewired at random with a probability p as in the 
one-dimensional model of Ref. [20]. The main differences with respect to the 
original model is that for any value of p we want to keep fixed the connectivity 
of each site. For such a reason we have implemented a rewiring procedure as in 
fig.2 in which the connections are rewired in e-"- 1 — w ~ „u„„„„ „ „;j.„ „• „„j 



Fig. 2. On the left, a schematic picture of the rewiring procedure to interpolate 
between a regular and a random topology by keeping fixed and equal to 4 the degree 
of each site. On the right, we report the characteristic path length L vs. the rewiring 
probability p. 

the edge %\ — % 2 that connects site i\ to its nearest neighbor % 2 in a clockwise 
sense. With probability p we decide whether to rewire this edge or to leave it 
in place. If the edge has to be rewired we (a) choose at random a second site 
ji and one of its edges, for instance the edge j\ — j 2 connecting site ji to site 
j 2 , and (b) we substitute the couple of edges i\ — i 2 and ji — j 2 with the couple 
ii -j 2 and j 1 - i 2 . 

We repeat this process by moving over the entire square lattice considering 
each site in turn until one lap is completed. In such a way the limit case p = 1 
is a QR graph with fixed connectivity (q) equal to 4. In the intermediate 
cases < p < 1 we can investigate the effects of an increasing number of 
long-range connections on the criticality of the model. Indeed, at a critical 
region of the parameter p between the regular (p = 0) and random (p = 1) 
networks, the topology produced by such a method exhibits a small-world 
behavior, characterized by the fact that the distance between any two sites 
on the graph is of the order of that for a random network and, at the same 
time, the concept of neighborhood is preserved, as for regular lattices. For this 
reason, we expect to obtain SOC in a small-world topology; the introduction 
of a few long-range edges create short-cuts that connect sites that otherwise 
would be much further apart. 




\ - 
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3.1 Earthquakes size distributions 



In our simulations the starting point for the construction of the SW network 
is a two-dimensional square lattice L x L with three different sizes: L = 32, 64 
and 128; the corresponding number of sites is N = L 2 . We have considered 
up to 10 9 earthquakes to obtain a good statistics for the earthquakes size dis- 
tribution Pn{s). In fig.3 we report the power law distributions resulting for 
N = 64 2 , a = 0.21 (non-conservative OFC model) and for two rappresentative 
values of the rewiring probability p (actually we made the simulations also for 
many other values of p in the range [0,0.1]). In the same figure we report also 
the comparison with the earthquakes size distribution for the dissipative OFC 
model on a scale- free network (that will be discussed in section 4). 




Iog10(s) 

Fig. 3. Earthquakes size distributions for the non-conservative OFC model (with 
a = 0.21) on rewired 64x64 lattice with two rappresentative values of rewiring 
probability p in the range [0,0.1]. A stretched-exponential function is used to fit the 
pdf's cutoffs. In the inset we plot the two exponents r and a as a function of the 
rewiring probability p. Here we report also the pdf for SF OFC model (see sec. 4). 

All the curves can be fitted by a stretched-exponential function P/v(s) = 
As~ T e~( s l^" , where s is the size of earthquakes, £ is the characteristic length 
and r and a are two exponents. We notice that, increasing more and more 
the rewiring probability, the power law is practically lost. This can be better 
exploited by plotting the value of the two exponents r and o as a function of p 
in the inset in Fig.3. Indeed one can expect stretching-exponential in various 
cases of stochastic processes where many length scales appear. We note also 
that, above the value p ~ 0.01 for which a suddenly approaches 1, the power 
law for the pdfs progressively disappears. In the next subsection we will show 
that the cut-off in the earthquakes probability distribution scales with the 
system size (the so called finite size scaling ansatz) only around this rewiring 
threshold. 
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3.2 Finite Size Scaling 



In order to characterize the critical behavior of the dissipative SW OFC model, 
a finite size scaling (FSS) ansatz is applied, i.e. P/v(s) ~ N' 13 f(s/N D ) where 
/ is a suitable scaling function and f3 and D are critical exponents describing 
the scaling of the distribution function. In Fig. 4 we consider a = 0.21 and 
a rewiring probability p ~ 0.006. We show the collapse of P/v(s) for three 
different values of N, namely N = 32 2 ,64 2 , 128 2 . The distribution Pn(s) sat- 
isfies the FSS hypothesis reasonably well, with universal critical coefficients 
with small rewiring probability, but, increasing p, as shown in the previous 
subsection, there is no power law at all. The critical exponent derived from 
the fit of Fig. 4 are (3 ~ 3.6 and D = 2. This result is in agreement with the 
FSS hypothesis implying that, for asymptotically large N, P/v(s) ~ s~ T with 
T = p/D ~ 1.8. 




Fig. 4. Finite Size Scaling for dissipative OFC model (with a = 0.21) on a 
small- world topology for three different values of N, namely N = 32 2 , 64 2 , 128 2 . 
The critical exponent derived from the fit are D = 2 and [3 ~ 3.6 and the rewiring 
probability is equal to 0.00586. 



Therefore, showing both power law behavior and FSS, the dissipative SW 
OFC model (in a restricted range of rewiring probability) clearly exhibits self- 
organized criticality (see Table 1). Let us point out that on the SW rewired 
topology the system behaves as in the compact square lattice, but the oc- 
currence of a small amount of long-range links disseminate the earthquakes 
over the network and the biggest earthquake size scales with the lattice. On 
the other hand, if we increase the amount of long-range links above a certain 
threshold (p ~ 0.006), the mechanism of synchronization is corrupted and the 
scaling behavior disappears. 
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Table 1 

In this table we list the SOC properties for OFC models on different topologies: a two 
dimensional nearest neighbor lattice (2D NN lattice), an annealed random neighbor 
(ARN) graph, a quenched random (QR) graph, a quenched random (QR+2) graph 
with two sites with coordination 3, a small-world (SW) network and a scale-free 
(SF) network. We always consider these models in the case a = 0.21 (dissipative 
regime) and with open boundary conditions. In particular, we report when there is 
power law and finite size scaling in earthquakes size distribution, according to each 
kind of topology. 

4 The OFC model on a scale-free network 



Finally we investigate criticality of the non- conservative OFC model on a scale- 
free network. It is an example of network displaying a small characteristic path 
length and a power-law distribution p(k) ~ A; -7 in the node connectivity k (de- 
gree). By using the preferential attachment growing procedure introduced by 
Barabasi and Albert [11], we start from m + 1 all to all connected nodes and 
at each time step we add a new node with m links. These m links point to old 
nodes with probability pi = ^ , where is the degree of the node i. This 

procedure allows a selection of the 7 exponent of the power law scaling in the 
degree distribution with 7 = 3 in the thermodynamic limit (N — > 00). Here 
we consider a scale-free network with 7 = 3 and A" = 1000. 
In this case, the toppling rule in Eqs. 1 must be modified to take into ac- 
count that different sites have a different coordination number qi. Each site 
consequently has a different which we determined by requiring that the 
total fraction a of the force transferred from the unstable site to the nearest- 
neighbor sites is constant in the system, i.e., ctj = a/qf, here we consider the 
case a = 0.21. 

We have found that there is no criticality in the system since there is no 
power law in the earthquakes size distribution, as shown in fig. 3. As previ- 
ously observed and in agreement with previous investigations [9,31,32], this 
result indicates that if the disorder is too strong the critical signatures disap- 
pear and the SOC behavior is destroyed (see Table 1). 
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5 Conclusions 



In conclusion, in this paper we have investigated the dissipative OFC model 
on small-world and scale-free networks. 

We have shown that, at variance with OFC models on other topologies which 
are critical only in the conservative case, the dissipative small-world OFC 
model clearly reaches a critical state characterized by power law behavior of 
earthquakes size distribution with finite size scaling of cut-offs. Indeed, in a 
lattice with a small number of rewired links the underlying spatial structure al- 
lows partial synchronization of distant blocks of the system. We think that this 
process could reproduce the long-range earthquakes dynamical correlations in 
the earth crust, according to the most recent geophysical observations. 
On the other hand, on a scale-free topology we do not observe SOC proper- 
ties. We expected this behavior because the connectivity is not fixed; so the 
dynamics is not synchronized, the disorder is too strong and the critical state 
is destroyed. As future directions, it seems interesting to better investigate the 
influence of topology and the role of disorder on the self-organized criticality 
properties of the OFC models. 
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